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Abstract 

We examine dark quantum superposition states of weakly bound Rb2 Feshbach 
molecules and tightly bound triplet Rbi molecules in the rovibrational ground state, 
created by subjecting a pure sample of Feshbach molecules in an optical lattice to 
a bichromatic Raman laser field. We analyze both experimentally and theoretically 
the creation and dynamics of these dark states. Coherent wavepacket oscillations 
of deeply bound molecules in lattice sites, as previously observed in Ref. [1], are 
suppressed due to laser-induced phase locking of molecular levels. This can be un- 
derstood as the appearance of a novel multilevel dark state. In addition, the exper- 
imental methods developed help to determine important properties of our coupled 
atom / laser system. 

1 Introduction 

Very recently, several groups have produced dense, ultracold ensembles of molecules that 
are deeply bound [HI IH [3l |l|5l and in a ro-vibrational ground state ID l2l 111 13 . This 
was achieved by binary association of alkali atoms in ultracold ensembles via two differ- 
ent pathways: (1) photoassociation |l6l|7]| and (2) magneto-association at Feshbach reso- 
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nances [|7l[8l combined with stimulated Raman adiabatic passage (STIRAP) BU, a special 
coherent optical transfer method. In contrast to photoassociation, magneto-association 
only produces weakly -bound Feshbach molecules [|7l[8l. STIRAP can then be used to 
transfer these weakly-bound molecules to the rovibrational ground state. This method is 
coherent, efficient, fast, reversible, and highly selective. STIRAP is based on a counter- 
intuitive light pulse sequence giving rise to a dynamically changing dark superposition 
state (Fig. [Di) 

\DS) = {a2\f)-^i\g))/^aj + Q.i (1) 

In this paper, we deliberately replace the efficient but complex STIRAP transfer of 
Ref. [Ill with a simple square laser pulse scheme. This reveals interesting fundamen- 
tal processes and dynamics in the coupled atom / laser system, that would otherwise be 
hidden. In addition, this procedure allows us to determine important properties and pa- 
rameters of our system and to check for consistency with our theoretical model. We study 
the creation and lifetime of dark superposition states that contain a sizeable fraction of 
deeply bound molecules. These molecules are held in a 3D optical lattice. Because the 
lattice potential is much shallower for the deeply bound molecules than for the Fesh- 
bach molecules, and because the transfer is fast, the deeply bound molecules coherently 
populate several Bloch bands. In contrast to Ref. [[ll, where similar circumstances lead to 
coherent oscillations in the lattice, oscillations are suppressed in the experiment described 
here due to phase locking of all quantum levels involved. A novel dark state appears which 
is a superposition of up to 8 quantum levels. We investigate the limiting conditions under 
which oscillations set in. 
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2 Experimental setup and initial preparation of molecules 

We carry out our dark state experiments with a 50 /im-size pure ensemble of 3 x 10^ 
weakly bound Rb2 Feshbach molecules. The molecules are trapped in the lowest Bloch 
band of a cubic 3D optical lattice with no more than a single molecule per lattice site [fTOll 
and an effective lattice filling factor of about 0.3. The lattice depth for the Feshbach 
molecules is 60 Er, where Er = n^fP' /Ima^ is the recoil energy, with m the mass of the 
molecules and a = 415.22 nm the lattice period. Such deep lattices suppress tunneling 
between different sites. A pure ensemble of Feshbach molecules has been produced as 
follows. We prepare a cold cloud of 6 x 10^ ^^Rb atoms that are either Bose condensed 
or nearlxi Bose condensed in a loffe-type magnetic trap with trap frequencies ft)x,j,z = 
271 X (7, 19,20) Hz. Within 100ms we adiabatically load the atoms into the 3D optical 
lattice. After turning off the magnetic trap, we flip the spins of our atoms from their initial 
state \F = l,m/7 = — 1) to \F = l,m^ = +l)by suddenly reversing the bias magnetic field 
of a few G. This spin state features a 210 mG-wide Feshbach resonance at 1007.40 G [[TT|. 
By adiabatically ramping over this resonance, we efficiently convert atoms at multiply 
occupied lattice sites into Rb2 Feshbach molecules. After conversion, inelastic collisions 
occur at lattice sites that contain more particles than a single Feshbach molecule, leading 
to vibrational relaxation of these molecules, release of binding energy into kinetic energy 
and removal of all particles from these sites. A subsequent combined microwave and 
optical purification pulse removes all remaining chemically unbound atoms, creating a 
pure sample of 3 x 10^ Feshbach molecules. Afterwards, the magnetic field is set to 
1005.8 G, where the Feshbach molecules are in a quantum state |/) which correlates with 
\F = 2, m/7 = 2, /i = 2, /2 = 2, V = 36, / = 0) at G. Here, F and /i^2 are the total angular 
momentum quantum numbers for the molecule and its atomic constituents, respectively. 



It turns out that this increases the number of Feshbach molecules. 




Internuclear distance (Bohr radii) 5/2n (IVIHz) 

Figure 1: (a) A- type three-level scheme for dark states. The lasers 1 and 2 couple the 
molecular levels |/), \g) to the excited level \e) with Rabi frequencies ^1^2, respectively, 
(b) Typical dark resonance. Shown is the fraction of Feshbach molecules |/) remaining 
after exposing them to both lasers in a 3 /is square pulse for varying two-photon detuning 
5. The data is identical to that shown in Ref. [JJ. 

and nif is the total magnetic quantum number; v is the vibrational quantum number for 
the triplet ground state potential (fl^S^ ) and I is the quantum number for rotation. 

The bichromatic Raman laser field for the creation of the molecular dark states is 
based on two lasers (1 and 2) which connect the Feshbach molecule level |/), via an 
excited level to the absolute lowest level in the triplet potential |^) (Fig.[T^). Laser 1 
is a Ti:Sapphire laser and laser 2 is a grating- stabilized diode laser. Both lasers are Pound- 
Drever-Hall locked to a single cavity which itself is locked to an atomic ^^Rb-line. From 
the lock error signals, we estimate frequency stabilities on a ms-timescale of 40 kHz and 
80 kHz for lasers 1 and 2, respectively. Both laser beams have a waist of 130 /im at the 
location of the molecular sample, propagate coUinearly, and are polarized parallel to the 
direction of the magnetic bias field. Thus, the lasers can only induce n transitions. 

The ground state |^) has a binding energy of 7.0383(2) THzx/z and can be described 
by the quantum numbers \F = 2, nip = 2,S = 1,7 = 3, v = 0,/ = 0) where S and / are 
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the total electronic and nuclear spins of the molecule, respectively. At 1005.8 G |^) is 
separated by hundreds of MHz from any other bound level, so that there is no ambiguity 
as to which level is addressed. The level \e) is located in the vibrational v = 13 manifold of 
the electronically excited ^Eg (5s + 5p) potential and has Ig character. It has an excitation 
energy of 294.6264(2) THz x h with respect to |/) , and a width T = 2;r x 8 MHz. The Rabi 
frequencies 2 of the two lasers depend on their respective intensities /i 2, i-C-, ^^i = 
2;r X 0.4 MHz a//i / (Wcm-2) and ^^2 = 2;rx30MHz ^Jh/i^cmr^), and are typically 
chosen to be in the MHz regime. 

3 Dark state evolution within a square pulse 

Our square pulse projection experiments are carried out as follows. We expose the Fes- 
hbach molecules |/) in the lattice to square pulses of Raman lasers 1 and 2 of variable 
pulse duration. Laser 2 is switched on about 1 /is before laser 1 to avoid excitation from 
I/) to \e) due to jitter in the laser pulse timing. The Raman lasers are resonant (5 = 0) 
and the Rabi frequency ^^2 ~ 2;r x 7 MHz while Q.\ is varied (Fig.[2l). After the pulse, we 
measure the fraction of molecules remaining in state |/) by dissociating them into pairs 
of atoms at the Feshbach resonance, releasing them from the lattice and applying standard 
absorption imaging. It is important to note that we actually only count atoms in the lowest 
Bloch band of the lattice. The release from the optical lattice is done as described in [[T2|. 
where after 13 ms of ballistic expansion we map out the Bloch bands in momentum space 
(see Appendix for details). 

Figure [2] shows the remaining fraction of molecules in state |/) versus pulse duration. 
Within 1 /is we observe a rapid loss of molecules that depends on the ratio ^2/^1- The 
remaining molecules are stable on a much longer timescale. This can be understood in 
terms of formation of a dark state \DS) . We can write 
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Figure 2: Dark state formation and lifetime. Shown is the fraction of Feshbach molecules 
remaining after subjecting them to a square pulse of Raman laser light of varying length 
for various Rabi frequency ratios ^^2/^1 (^^2 ~ 2;r x 7 MHz). After switching on the 
lasers, a certain fraction of molecules is lost within 1 /is and a dark state has formed 
which has a much longer lifetime. The solid lines represent model calculations (Sec. Sj) 
which can be used to determine the Rabi frequencies and short-term laser linewidths. 



I/) = {Q2\DS) + ^ii\BS))/^al + al (2) 

where 

\BS) = {ai\f)+a2\g))/^^]+nl (3) 

is a bright state which quickly decays via resonant excitation to level \e). The dark 
state remains after the lasers are switched on and can be detected as a fraction Q.i/(Q.^ + 
0.2)'^ of molecules projected back to |/) after switching off the lasero 

Also, after the pulse a fraction Q.\0.2/ + ^2)^ of the initial molecules are in state 
\g) with a maximum of 25% for Q.i = Q.2. Thus, a sizeable fraction of the molecules can 

■^This fact can be used to conveniently calibrate the Rabi frequency ratio Q.i/Q.2. We found good agree- 
ment with other calibration methods for the Rabi frequencies. 
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be coherently transferred to the ground state. Remarkably, this transfer takes place in less 
than 1 /is! Such short transfer times cause Fourier broadening, resulting in considerably 
reduced laser stability requirements. In addition, due to the formation of a dark state, 
there is still a well-defined phase relation between the |/) and \g) molecules. 

As can be seen from Fig. [2l the dark state slowly decays. Its lifetime is shortest for 

= Q.2, where we measure it to be ~50 /is. The decay of the dark state is likely due 
to phase fluctuations of the Raman lasers. Phase fluctuations lead to an admixture of a 
bright state component to the otherwise dark state, which causes losses. In Sec.|4]we will 
show that these fluctuations can be expressed in terms of the short-term relative linewidth 
of the lasers, 7, which we find to be about In x 20 kHz. In principle, the decay of the dark 
state could be due to other effects, such as coupling to levels other than |/), \e), and |^). 
However, we have verified that this is not the case, because losses due to optical excitation 
are completely negligible on the 100 /is-timescale when we expose a pure ensemble of |/) 
(1^)) molecules to only laser 2(1). 

We also searched for laser power dependent shifts of the two-photon resonance. Using 
the Raman square pulse measurements, we scanned the relative detuning of the lasers for 
a fixed pulse duration and various laser powers. Within the accuracy of our measurements 
of In X 200 kHz, we could not detect any shifts of the resonance. 

The behavior in Fig.[2]is described well by a closed three-level model (a A system) and 
its dynamics can be simulated with a master equation which we describe in the following. 

4 Three-level model and master equation 

Neglecting lattice effects, we can describe the internal dynamics of the molecules as they 
are subjected to the Raman laser fields with a three-level model. We use a master equa- 
tion [dSldll which takes into account decoherence due to phase fluctuations of the Raman 
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Figure 3: Level scheme for the master equation. 

lasers. We consider the case where laser 1 is kept on resonance and laser 2 has a detun- 
ing d (Fig. |3]). Identifying the levels |/), \g), \e) with numbers 1, 2, 3, respectively, we 
can write the master equation as, 



dp 
dt 



.3^ 



3yt 



k=\ 



_ V (^33.^33. p^p. ^33. ^33) 



(4) 



where p is the density matrix, Qi 2 are the Rabi frequencies, Y is the spontaneous 
decay rate of the excited level and 7 is the relative linewidth of the two Raman lasers. 
The matrices cr" and cr" are ladder operators and each is the transpose of the other. For 
example 



.32 




(5) 



Setting the linewidth of the excited level F = 8 MHz, the detuning 5 = and Rabi 
frequencies = 2;r x 7 MHz and ^\ to give the ratios in Fig.[2l we fit all the data with 
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Figure 4: We plot the transfer efficiency for the round-trip STIRAP process as a function 
of the hold time T/, between the two STIRAP pulses. With our procedure we only count 
molecules whose constituent atoms end up in the lowest Bloch band after transfer. The 
oscillations in the transfer efficiency are due to breathing oscillations of localized spatial 
wavepackets of molecules in the lattice sites. The solid line is from a multi-band model 
calculation (Sec.[6l). This plot is taken from Ref. 

a single fit parameter 7. As a best fit, we obtain a relative linewidth of the two Raman 
lasers 7 = 2;r x 20 kHz, which is a reasonable value for our laser system. 



5 Coherent oscillations and their suppression 

In reference |[T1 coherent oscillations of molecular wavepackets of |^) molecules in the 
optical lattice were observed. We now investigate how these observations fit together 
with the experimental results of the square pulse projection experiments presented here. 
For clarity, the oscillation data from Ref. [IJ are presented again in Fig. |4]and briefly 
discussed. 

Using an STIRAP pulse sequence, Feshbach molecules are efficiently transferred to 
level 1^). The Raman lasers are extinguished and the molecules are held for a time T/„ 
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Figure 5: Wavepacket dynamics. Directly after the STIRAP transfer of a molecule from 
I/) to 1^) the shape of its wavepacket (thick solid line a) essentially corresponds to the 
vibrational ground state of the sinusoidal lattice potential for Feshbach molecules Vf (thin 
solid line). In the much weaker potential felt by the ground state molecules Vg (thin dashed 
line) the wavepacket starts to oscillate. After 1/8 of the oscillation period Xosc = In I Oh 
its shape roughly corresponds to the vibrational ground state for |^) (thick dashed line j8) 
and reaches its maximum extension after XosJ^ (thick dash-dotted line 7). 

after which they are transferred back to |/) with a reverse STIRAP pulse. The number 
of recovered Feshbach molecules is counted. However, we only detect atoms that end up 
in the lowest Bloch band after dissociation of the Feshbach molecules (see Appendix). 
The oscillation can be understood as follows. We consider the localized spatial center- 
of-mass (c.o.m.) wavepacket of a Feshbach molecule at a particular lattice site in the 
lowest Bloch band. The first STIRAP transfer projects this wavepacket onto the much 
shallower lattice potential felt by the \g) molecules (Fig. [5]) without changing its shape. 
As a consequence, |^) molecules are coherently spread over various Bloch bands, and 
the wavepacket undergoes "breathing" oscillations with the lattice site trap frequency w,. 



^Due to a smaller dynamic polarizability, the lattice depth for the tightly bound |g) molecules is shal- 
lower than for the Feshbach molecules by a factor of wlO. 
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These coherent oscillations (period ^SO fis ) are damped by tunneling of |^) molecules 
in higher Bloch bands to neighboring lattice sites. The reverse STIRAP transfer maps 
this periodic oscillation back to the Feshbach molecule signal in Fig. IH Higher Bloch 
bands are populated here as well, but are at most partially counted in our scheme (see 
Appendix), which leads to an apparent decrease in transfer efficiency. 

The question arises why similar oscillations are not observed in our square pulse pro- 
jection measurements shown in Fig. [21 especially for the case Q.i = Q2 where 50% of 
the population is in state \g). One might assume that the spatial wavepackets of the |^) 
molecules undergo similar breathing oscillations. These oscillations would then periodi- 
cally break up the dark superposition state and lead to corresponding losses. They would 
also periodically produce population in higher Bloch bands of the Feshbach molecule lat- 
tice. As we will see, the oscillations are suppressed because the Raman lasers phase lock 
the involved quantum levels which stops, in a sense, the free evolution of the wavepack- 
ets. We can understand this behavior in detail with the help of a multi-band model, which 
we describe in the following. 

6 Multi-band model 

In an optical lattice the molecular levels |/), \g) and \e) from the previous model have 
a substructure given by the lattice Bloch bands. Because the lattice depths for the levels 
I/), 1^) and \e) are in general different, the respective band structures will also vary. This 
combination of external (c.o.m. motion in the lattice) and internal degrees of freedom 
gives rise to a number of new quantum levels which are coupled by the laser fields (Fig. [6]). 
We assume each Feshbach molecule to be initially localized in a singly-occupied lattice 
site. The corresponding localized molecular wavepacket can be described by Wannier 
functions [,15 J which form a complete set of orthonormal functions. In the following we 
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Figure 6: Multi-band model. The three molecular levels |/), \e) and |^) have a Bloch 
band substructure due to the optical lattice. We restrict the model to the 4 lowest Bloch 
bands with even symmetry (band index n = 0, 2, 4, 6). 

will denote the Wannier function for level | a) and band n as \^an)- We note that for deep 
lattices, these Wannier functions closely resemble harmonic oscillator wavefunctions. 

The Raman lasers couple different |^a„) according to the respective wavefunction 
overlaps (Fig. (6]). Since the initial wavepackets of the Feshbach molecules are symmetric, 
only even bands will be populated. We restrict our calculations to the four lowest Bloch 
bands with even symmetryjj corresponding to the band indices n = 0, 2, 4, 6. The dynam- 
ics in each of the three lattice directions is then described by a 12-level model, which can 
in principle be solved in terms of a master equation (Sec. Sj). However, we have used a 
Schrodinger equation-based model since the numerical code is less involved. In this ap- 
proach, laser phase fluctuations are neglected, and we introduce a lattice site tunnel rate 
for each band. These tunnel rates are chosen to match the expected tunnel rates for the 
different bands and are slightly adjusted for a better fit of the data in Fig. HI We note that 
the results of the model calculations are essentially independent of the excited state lattice 
depth, which is not well known. 



"^The effect of including higher bands into the model was found to be negligible. 
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The Hamiltonian H of our time dependent Schrodinger equation 

in m = Hm (6) 

at 

has the form of a 12 x 12 matrix, 

/ Efo-^Jfo ^ai{t)-Mfo^,o ...\ 

H = h j^l{t)-MeOJO ^a2{t) ■ MeO,gO EeQ-^T-^M \^l{t)-MMJ2 ••• 
\^i{t)-Mf2,eQ Ef2-{jf2 ... 

\ \ \ \ \ I 

(7) 

Here E^n and Jan are the energy and tunnel matrix element respectively for the Wan- 
nier function |^a„) in band n of level M^^y^^^y- = {^'anl^pk) the overlap integral of 
the respective Wannier functions. 

Diagonalizing this Hamiltonian, we find twelve "eigenstates" of the coupled system 
which in general have complex eigenvalues. In the following, we study the case of strong 
coupling (^2i 2 ^ w?) y which is the regime for phase locking. In this regime, four of these 
eigenstates have negligible contribution from the exited level \e) and thus a long lifetime. 
These 4 quasi-dark states essentially correspond to the 4 lattice bands in our model and 
will be denoted as \DSn) with n = 0, 2, 4, 6. We now study the spatial waveforms of these 
dark states (Fig.|7]) and compare the components with |/) and |^) character. Neglecting a 
small \e) component the dark superposition state, \DSn) can be written as 

\DSn) = \g){g\DSn) + \f){f\DS„). (8) 

As an example = Q.2) Fig. |7] shows that the wavepackets {g\DSn) and {f\DSn) have 
the same shape. This is not surprising since this ensures that the ratio of the |/) and |^) 
amplitudes equals ^22/^1 everywhere, as in Eq.[T] 



For our experiments where Q.12 ^ 27r x 1 MHz and (Ot ^27t x lOkHz this condition is satisfied. 
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Figure 7: Spatial wavepackets of the dark state \DSo) and the Wannier functions \^fo), 
l^^o)- The dark state \DSo) has two components, one having |/) character {{f\DSn)) and 
the other one having |^) character {{g\DSn)). The wavepackets of and {g\DSn) 

essentially have the same shape. They are mainly composed of the lattice ground states 
\^fo) (thin solid line) and |^go) (thin dotted line). All depicted states are normalized. 
The parameters used are = ^2 = 27r x 7 MHz, Vf = 60 Er and Vg = 6Er, as in our 
experiments. 

Let us now discuss the formation and evolution of the dark state that we have observed 
in the square pulse experiments of Sec.[3l A dark state \DS) is formed in less than 1 fls by 
subjecting Feshbach molecules to a square Raman laser pulse. As in the STIRAP transfer 
(discussed in Sec. |5]) the initial projection onto \DS) will not change the shape of the 
Feshbach molecule wavepacket, given by the Wannier function |^/o)- The dark state can 
be expressed as a coherent superposition of the four dark eigenstates \DSn} of the 12-level 
Hamiltonian 



The subsequent coherent evolution of these dark states will again in principle lead to 
breathing oscillations. The amplitude of these oscillations depends on the extent to which 
higher bands (i. e., \DS)n, n > 0) are excited. The excitation increases with increasing 
deviation of \DS) from the initial state |/), i. e., with rising Q.1/Q.2. 



\DS)= I 



Cn\DSn). 



(9) 



n=0,2,4,6 
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This can also be understood from another point of view. The effective lattice potential 
felt by the molecules in such a superposition state is the weighted average of the potentials 
for the two contributing states |/) and |^). For the case Q.i = Q.2 this effective potential 
is about half as deep as that for the Feshbach molecules. Compared to the case of pure 
ground state molecules (Fig. H]) where the lattice potential is reduced by a factor of 10, 
the oscillations of the wavepacket are strongly suppressed and cannot be observed with 
our current experimental precision. For Q.i ^ Q.2, the dark state \DS) has a dominant 
contribution from state |^), and the effective lattice potential essentially corresponds to 
the one for ground state molecules. In this case oscillations appear despite the strong 
coupling, a fact which we also have experimentally verifieq^. 

Conclusion We have analyzed coherent wavepacket dynamics and their suppression in 
a 3D optical lattice. We observed optically induced phase locking of a number of quantum 
levels, which can also be viewed as the appearance of a novel multi-level dark state. The 
experiments were carried out with tightly bound molecules as a component of a dark 
quantum superposition state. Thus, the experiments demonstrate control of molecular 
motion in an optical lattice for the first time. In addition, different models have been 
introduced and discussed in detail, with which the lattice dynamics can be understood 
and quantitatively described. 
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support in model calculations. We also thank Gregor Thalhammer for early assistance 
in the lab, and Florian Schreck for loaning us a Verdi VI 8 pump laser. This work was 
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^For these experiments we have ramped into the dark state and back in a fashion similar to STIRAP 
transfer pulses to avoid strong losses caused by direct projection into \DS) w |^) {g\DS). 
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Figure 8: Shown is a typical absorption image which displays the atomic quasi- 
momentum distribution in the optical lattice after exposure to the Raman laser beams 
and subsequent adiabatic molecule dissociation. Atoms inside the square region come 
from the lowest Bloch band. Infi/a is the modulus of the reciprocal lattice vector. 

Appendix: Theoretical band population analysis 

As stated before, our signals only include molecules for which the constituent atoms end 
up in the lowest Bloch band of the lattice. A controlled lattice rampdown in a few mil- 
liseconds maps the bands and quasi-momentum distribution of the atoms into momentum 
space [fT6l[T2l . We image these distributions after 13 ms of time- of- flight via absorption 
imaging. Fig. [8] shows a typical distribution. The dotted square region corresponds to the 
lowest Bloch band and is dominantly populated. 

An important question is how the Bloch bands for the Feshbach molecules map onto 
the Bloch bands for the atoms. In other words, if we measure the atomic population of the 
Bloch bands - do we know what the band population for the molecules was? As the lattice 
is very deep for the Feshbach molecules and atoms, we can approximate the potential 
at an individual lattice site as harmonic with trap frequency ft)/. In one dimension, the 
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eigenfunctions of the harmonic oscillator are 



1 / I f x\^\ f X 



\^n) = , , ^ exp -- - HJ - , (10) 

where xq = ^Jfij (Otm is the oscillator length and is the n* Hermite polynomial. 
We assume that we have two atoms in a lattice site with coordinates xi 2- The relative and 
c.o.m. coordinates of the atom pair are 



Xr = 1/V2(JC1-JC2) (11) 

Xc = \l\fl{x\^xj) (12) 

The c.o.m. (U(. In one dimension, the eigenfunctions of the harmonic oscillator are We 
assume that we have two atoms in a lattice site with coordinates xi 2- The relative and 
c.o.m. coordinates of the atom pair are 



Xr = \ I y 7.{x\ — xt) (13) 
Xc = 1/V2(jci+;C2) (14) 

The c.o.m. potential Vc for the pair will be harmonic with trap frequency cOi and the 
potential Vr for the relative coordinate will be a sum of the harmonic potential and the 
interaction potential (Fig. ID). 

When we form or dissociate a molecule by adiabatically ramping across a Feshbach 
resonance, only the quantum level in the Vr potential will change - from a molecular 
bound state to an unbound atomic pair state in the lowest Bloch band. The wavefunction 
in the c.o.m. coordinate remains unchanged. We can now calculate how band populations 
of Feshbach molecules converted to atomic band populations by using the coordinate 
transformations. As an example: A Feshbach molecule in the lowest Bloch band (i. e., 
center-of-mass coordinate) will produce an atom pair with the following wavefunction: 
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Figure 9: Potentials for the center-of-mass and relative coordinate of two atoms trapped 
at a site of the optical lattice, (a) The bound states (solid lines) of the harmonic center- 
of-mass potential Vc correspond to the molecular Bloch bands, (b) At short interatomic 
distances the relative potential Vr is dominated by the interaction potential, which allows 
the formation of bound molecular states (solid lines). Unbound atoms are trapped by the 
lattice potential at larger separation (dashed lines). With the help of a Feshbach reso- 
nance the lowest trap state can be converted into a high molecular state. Note that in this 
schematic view both energy and distance for the two contributions to the relative potential 
are not to scale. 

1^) °^ exp(-l/2 x^) exp(-l/2 xj) = exp(-l/2 xf) exp(-l/2 xj). This means that both 
atoms will also end up in the lowest Bloch band of the lattice. This analysis can be ex- 
tended to any band. Table \T\ gives the conversion amplitudes from molecular to atomic 
bands for the four lowest symmetric molecular bands. Correlations between the two con- 
stituent atoms of a molecule are not discussed here. 

We finally note that when we apply absorption imaging, the optical density of the 
atomic sample is integrated in the direction of observation. Thus in this direction no band 
population analysis is possible. We accounted for this in our multi-band model described 



in Sec. |6l 
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Table 1 : Band conversion amplitudes in the harmonic oscillator approximation. Each line 
gives the amplitudes for a constituent atom of a molecule in a certain band to populate 
various atomic bands after dissociation. 
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